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FUgy mfiic ftiiinpn iMvgRjsg trpatmfmt PLAMMiMra 
BACKGROUND OF THE INVENTION 

10 TECHNICAL FIELD 

The present invention relates to inverse treatment planning. More specifically, 
the present invention relates to an artificial intelligence method for guiding inverse 
treatment planning. 

15 BACKGROUND ART 

Modern day radiation therapy of tumors has two goals: eradication of the tumor 
and avoidance of damage to healthy tissue and organs present near the tumor. It is 
known that a vast majority of tumors can be eradicated completely if a sufficient 
radiation dose is delivered to the tumor; however, complications may result from use of 
20 the necessary effective radiation dose, due to damage to healthy tissue which 
surrounds the tumor, or to other healthy body organs located close to the tumor. The 
goal of conformal radiation therapy is to confine the delivered radiation dose to only the 
tumor volume defined by the outer surfaces of the tumor, while minimizing the dose of 
radiation applied to surrounding healthy tissue or adjacent healthy organs. 

25 

Conformal radiation therapy has been traditionally approached through a range 
of techniques, and typically uses a linear accelerator ("LINAC") as the source of the 
radiation beam used to treat the tumor. The linear accelerator typically has a radiation 
beam source that is rotated about the patient and directs the radiation beam toward the 
30 tumor to be treated. The beam intensity of the radiation beam is predetermined, 
constant beam intensity, ly^ultileaf collimators, which have multiple leaf, or finger. 



projections that can be moved individually into and out of the path of the radiation 
beam, can be programmed to follow the spatial contour of the tumor as seen by the 
radiation beam as it passes through the tumor, or the °beam's eye view* of the tumor 
during the rotation of the radiation beam source, which Is mounted on a rotatable gantry 
5 of the linear accelerator. The multiple leaves of the multileaf collimator form an outline 
of the tumor shape as presented by the tumor volume in the direction of the path of 
travel of the radiation beam, and thus block the transmission of radiation to tissue 
disposed outside the tumor's spatial outline as presented to the radiation beam, 
dependent upon the beam's particular radial orientation with respect to the tumor 
10 volume. 

Another approach to conformal radiation therapy involves the use of 
independently controlled collimator jaws that can scan a slit field across a stationary 
patient at the same time that a separate set of collimator jaws follows the target volume 
15 as the gantry of the linear accelerator rotates. An additional approach has been the use 
of attachments for LINACs that allow a slit to be scanned across the patient, the 
intensity of the radiation beam in the entire slit being modified as the slit is being 
scanned. 

20 A further approach for conformal radiation therapy treatment has been the use of 

a narrow pencil beam of high energy photons, whose energy can be varied, and the 
beam Is scanned over the tumor target volume so as to deliver the best possible 
radiation dose distribution in each orientation of the gantry upon which the photon beam 
source is mounted. 

25 

A major problem associated with such prior art methods of conformal radiation 
therapy are that if the tumor volume has concave borders, or surfaces, varying the 
spatial configuration, or contour, of the radiation beam, the therapy is only successful 
part of the time. In particular, when the convolutions, or outer surfaces, of a tumor are 
30 re-entrant, or concave, in a plane parallel to the path of the radiation treatment beam, 
healthy tissue or organs may be disposed within the concavities formed by the outer 
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tumor cx)ncave surfaces, as well as the fact that the thickness of the tumor varies along 
the path of the radiation beam. 

In order to be able to treat tumors having concave borders, it is necessary to vary 
5 the intensity of the radiation beam across the surface of the tumor, as well as vary the 
outer configuration of the beam to conform to the shape of the tumor presented to the 
radiation beam. The beam intensity of each radiation beam segment should be able to 
be modulated to have a beam intensity related to the thickness of the portion of the 
tumor through which the radiation beam passes. For example, where the radiation 
10 beam is to pass through a thick section of a tumor, the beam intensity should be higher 
than when the radiation beam passes through a thin section of the tumor. 

Dedicated scanning beam therapy machines have been developed wherein 
beam intensity modulation can be accomplished through the use of a scanning pencil 

15 beam of high energy photons. The beam Intensity of this device is modulated by 
increasing the power of its electron gun generating the beam. The power increase is 
directed under computer control, as the gun is steered around the tumor by moving the 
gantry upon which it is mounted and the table upon which the patient lies. The effect is 
one of progressively "painting" the target with the thickness, or intensity, of the paint, or 

20 radiation beam intensity, being varied by the amount of paint on the brush, or how 
much power is applied to the electron gun, as the electron gun moves over the tumor. 
Such dedicated scanning beam therapy machines, which utilize direct beam energy 
modulation, are expensive and quite time consuming in their use and operation, and 
are believed to have associated with them a significant patient liability due to concerns 

25 over the computer control of the treatment beam itself. 

Other methods and apparatus for conformal radiation therapy have been 
developed that spatially modulate the beam intensity of a radiation beam across a 
volume of tissue in accordance with the thickness of the tumor in the volume of tissue 
30 by utilizing a plurality of radiation beam segments. Such methods and apparatus utilize 
attenuating leaves, or shutters, in a rack positioned within the radiation beam before the 
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beam enters the patient. The tumor is exposed to radiation in slices, each slice being 
selectively segmented by the shutters. However, a minor disadvantage of that method 
and apparatus results from the fact that only two slices of tissue volume may be treated 
with one rotation of the gantry of the linear accelerator. Mthoug^ the slices may be of 
5 arbitrary thickness, greater resolution is accomplished by selecting slices for treatment 
that are as thin as possible. As the thickness of the treatment slices decreases, the time 
it takes to treat the patient increases because more treatment slices are required in 
order to treat the entire tumor volume. 

10 The foregoing methods and apparatus are designed to minimize the portion of 

the structures being exposed to radiatk>n. However, because exposure to surrounding 
structures cannot be completely prevented, treatment plans are desired that are 
optimized to eradicate the tumor volume while minimizing the amounts of radiation 
delivered to the surrounding structures. Existing methods and apparatus for optimizing 

15 treatment plans use a computer to rate possible plans based on score functions which 
simulate a physician's assessment of a treatment plan. However, existing methods and 
apparatus have proven to be insufficient. 

Existing methods and apparatus utilize a computational method of establishing 
20 optimized treatment plans based on an objective cost function that attributes costs of 
radiation of various portions of both the tumor and sunrounding tissues, or structures. 
One such computational method is known in the art as simulated annealing. Existing 
simulated annealing methods utilize cost functions that consider the costs of under- 
exposure of tumor volumes relative to over-exposure of surrounding structures. 
25 However, the cost functions used in existing methods do not account for the structure 
volumes as a whole, relying merely on costs related to discrete points within the 
structure, and further do not account for the relative importance of varying surrounding 
structure types. For example, certain structure types are redundant in their function and 
substantial portions of the structure volume can be completely eradicated while 
3 0 retaining their function. Other structure types lose their function if any of the structure is 
completely eradicated. Therefore, the more sensitive structure volumes can receive a 
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measured dose of radiation so long as no portion of the structure is subjected to a letlial 
dose. 

Existing cost functions utilized in the optimization of treatment plans do not 
5 account for such varying costs associated with the different types of structures. After 
the treatment plan is optimized, the physician currently must evaluate each computed 
treatment plan for compliance with the desired treatment objective. If the computed 
treatment plan does not successfully meet the treatment objectives, the optimization 
process is repeated until a treatment plan can be computed that meets the physician's 
10 treatment objectives for both the tumor volume and the surrounding structures. Further, 
existing methods and apparatus do not allow the physician to utilize the familiar 
Cumulative Dose Volume Histogram ("COVH") cun/es in establishing the desired dose 
distributions. 

15 Prior to the development of the present invention, there has been no method or 

apparatus for conformal radiation therapy, for use with a radiation beam having a 
predetemilned, constant beam intensity for treatment of a tumor which: is simple and 
economical to use; that has what is believed to be a high safety factor for patient safety; 
computes an optimal treatment plan to meet conflicting, pre-detemiined, treatment 

20 objectives of a physician, accounting for objectives in both the target tumor volume and 
multiple structure types; and utilizes provides the desired dose distributions for each 
target tumor volume and tissue and structure types. 

DESCRIPTION OF THE INVENTION 

25 An artificial intelligence (A!) method, fuzzy logic, is applied to optimize 

parameters In the Inverse treatment planning for intensity-modulated radiation therapy 
(IMRT). With the capability of fuzzy inference, the parameter modification of the 
objective function is guided by physician's treatment intention and experience. For the 
different parameters involving inverse planning, the corresponding fuzzy Inference 

30 systems (FISs) are developed in order to accomplish the treatment requirement. With 
the use of fuzzy inference, the efficiency and quality of inverse planning can be 
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substantially improved. 



The present invention provides a method of using fuzzy logic to optimize 
treatments of patients. More specifically, the present invention uses a fuzzy inference 
5 systems that uses three modules: a Fuzzifier, a Inference engine, that includes fuzzy 
rules, and a Defuzzifier. During the process of fuzzification a single input value is 
compared to the membership functions as defined for that input variable. If the 
membership functions have a nonzero output it will take effect in the final results of the 
FIS. The Fuzzifier calculates the response of rules for the input values and the 
10 inference engine modifies the consequent rules in response to input values. The 
Defuzzifier generates a final output based on the result of the inference engine. 

Preferably, a computer is utilized to automatically calculate this in response to 
the input of data, however, a human can also calculate the same using the formula 
15 show in Appendix A. A detailed description of formula is set forth in Appendix A 
included herewith and incorporated by reference in its entirety. 
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IP application for 
artificial intelligence guided inverse planning system 

Title; Artificial Intelligence-guided Inverse Treatment Planning 

Abstract: An artificial intelligence (AI) method, fuzzy logic, is applied to optimize 
parameters in the inverse treatment planning for intensity-modulated radiation 
therapy (IMRT). With the capability of fuzzy inference, the parameter modification 
of the objective fiinction is guided by physician's treatment intention and experience. 
For the different parameters involving inverse planning, the corresponding flizzy 
inference systems (FISs) are developed m order to accomplish the treatment 
requirement. With the function of fiizzy inference, the efficiency and quality of 
inverse planning will be substantially improved. 

Introduction: 

The conventional inverse treatment planning involves optimizing an intensity 
spectrum by an objective function in order to obtain an ideal dose distribution. 
Sometimes, the dose distribution calculated by the optimized intensity spectrum might 
not be the one we expected. It is often resulted by an improper selection of parameters 
selected in the objective function. As the desired dose distribution can only be achieved 
by a refined combination of parameters that are usually unknown, it was inevitable for a 
planner to try different combinations of parameters in order to achieve a desired dose 
distribution. At present, this process is done by a trial-and-error approach.^^'^^^ It is not 
only time-consuming but also difficult to find a desired combination of parameters. 

A number of methods have been proposed in recent years in order to tackle this 
problem.^^*^'^^"^^' Starkschall^^^^ proposed an approach that removed the necessity of 
defming a *T)est" treatment plan, and incorporated the dose-volume constraints into a 
system to search for a feasible plan that could satisfy the constraints. If no calculated 
doses satisfy the treatment goal, the planner will provide a guide about how the dose- 
volume constraints may be modified to achieve a feasible result. This approach is only 
applied to die conventional three-dimensional (3D) treatment planning. Wu and Mohan^*' 
developed an optimization system, which employed both dose- and dose-volume-based 
objective functions. In their system, the optimal plan is selected by calculating the cost of 
the objective fimction, or "plan score" (the lower the score, the better the plan). Xing et 
al.^^^*^"^' presented a method which employed a second stage evaluation function to 
compute the differences between the calculated and the ideal dose volume histograms. 
Based on the results of the evaluation fimction, the weighting factors in the objective 
function are adjusted. This procedure is to minimize both the objective and evaluation 
functions in a round-robin manner. Later, further improvement is achieved by using a 
statistical measure called preference function, which is constructed based on the 
empirical judgment. The problem about the selection of weighting factor still exists 
because it translated to the problem of how to specify the parameters in the evaluation or 
preference function. A similar method was also proposed by Wu et al.^*^^ using a genetic 
algorithm to optimize the weighting factors and beam weights in the conventional 3D 



treatment planning. Li and Yin^^^ introduced fuzzy logic into the inverse planning system 
to adjust the weighting factors for normal tissue. The result was promising. However, 
optimizing the parameters for the target and critical organ were not included in the 
system. Also, the weighting factors initialized by the fiizzy functions still need to be 
modified by the trial-and-enor approach. 

Purpose: 

One of the major advantages of IMRT is to spare dose of the critical organ 
adjacent to the target. When a strict compromise between target dose and critical dose 
was required, fine-tuning of parameters in the objective function is important and 
necessary. As the effect of parameters on the output of the inverse planning is uncertain, 
it is hard to say which combination of parameters is best for the current case. Generally, 
this decision is made according the priority of organs involved and is individualized. 
More concerns will pay to the organ with high priority. As the trial-and-eiror approach 
can be regarded as a way using empirical knowledge to guide the parameter modification, 
it can be substitute by existed AI approaches, such as expert system and fuzzy inference 
system. 

The main purpose of this research is to use fuzzy interference system to defme 
parameters in the objective function based on treatment intention instead of trial-and- 
error approach. 

Detailed description 

See attached manuscripts. 

The advantage provided bv the invention 

Will improve the accuracy and efficiency of inverse treatment planning for IMRT. 
Specifically, 

(1) The heuristic and practical experience (firom physician, physicist, planner) can be 
used to guide the optimizing the parameters of inverse planning in order to improve 
the dose distribution; 

(2) The conformity of target dose distribution can be improved and high target dose 
could be feasible - improve the quality of inverse planning; 

(3) The time spent on trial-and-enor test can be reduced significantly and the planner will 
be free from this time-consuming task - improve the efficiency of inverse planning. 

Commercial significance of that advantage: 

It will make the treatment planning system more competitive. 

The closest known prior publications and prior uses, if anv. 
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SECnON 2 - RESEARCH PLAN 




A> SPECIFIC AIMS; 

The objective of this proposal is to develop an optimal approach for inverse treatment 
planning and a three-dimensional verification method for intensity-modulated radiation therapy. 
Specifically, this study will test two hypotheses: 

(1) The vague prescription could be represented by the theory of fuzzy set in inverse treatment 
planning. 

(2) The dose distribution delivered by limited number of static intensity-modulated radiation 
fields can be directly verified using a three-dimensional reconstruction method. 

The specific approaches to test these hypotheses are to: 

(1) Develop an effective inverse planning method to test prescription method and generate 
intensity spectrum for intensity-modulated radiation tiierapy. 

(2) Develop a logical fuzzy function to represent the vague dose prescription in the 
inverse-planning algorithm. 

(3) Develop a three-dimensional CT reconstruction algorithm for patient anatomy using a 
hmited number of projections. 

(4) Develop a method to reconstruct dose distribution in 3-D CT images based on a limited 
number of static IMRT fields using multilevel algebraic reconstruction technique. 

B. SIGNIFICANCE: 

It has been shown that the tumor loco-regional control can be improved with higher 
radiation dose than conventional prescription (Ref.?). However, delivery of higher dose is 
typically associated with the increase of probability of normal tissue complications. One approach 
to reduce this complication is to minimize the volume of normal tissue irradiated (ref.?). To 
achieve this goal, it will need a confonnal dose distribution to target volume, a dedicated dose 
delivery mechanism, and a verification method for patient positioning and dose delivered. 

Recent studies indicated that conformal dose distribution could be effectively achieved 
with the treatment technique called intensity-modulated radiation therapy (IMRT) (Ref?). Several 
promising delivery devices have also become available, such as static or dynamic MLC and 
tomotherapy to deliver conformal radiation dose (Ref. ?). The basic concept of IMRT is that a 
dedicated delivery device with an intensity-variable modulates a uniform intensity in a traditional 
treatment field. However, it is still a very challenging issue in terms of how to generate an 
effective and optimal intensity spectrum and how to verify modulated radiation delivery. The fust 
issue is also related to the problem of inverse treatment planning (or treatment planning 
optimization). 

Inverse planning method' describes a specific treatment planning procedur<e in which, 
differing from traditional approach, both dose and volume are given first. - Then a set of modulated 
beams is generated through a computer-aided optimization process in order to satisfy the 
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prescription. This process is extremely important if the shapes of the target arid critical organs are 
complicated, especially when the target has concavity and a critical organ lies in the hollow of the 
concavity. Typically, inverse treatment planning for intensity modulated radiation therapy 
involves the selection of an objective function and method of optimization. For a given objective 
function, an optimal treatment plan usually requires the optimization of beam intensity elements, a 
prescription method, and beam number and orientation. 

One of the most challenging problems in the optimization of treatment planning is how to 
construct a model by which the aim of radiation therapy can be fulfilled. The models that have 
been studied in the past can be classified as either physical or biological. There have been detailed 
discussions in recent literature (Ref 1,11) concerning the merits and lunitations of these two types 
of models. While biological models may be able to directly measure the clinical outcome, they still 
remain in the formative stages and suffer from controversy concerning the validity of the 
radiobiological response data used (such as, tumor control probability (TCP) and normal tissue 
complication probability (NTCP)) (Ref. 1, 13-15). On the other hand, the physical dosimetric 
prescription has been well established as the clinical norm. In the traditional physical models, one 
optimizes an objective function that is the measure of closeness of the calculated dose distribution 
to the prescribed dose distribution. The crucial problem here is how to give the optimal dose value 
for the normal tissue so that the two objectives, delivering the desired dose to the target volume 
and minimizing the dose to normal tissues, can be achieved accordingly. 

A quadratic model has typically being used in inverse treatment plaiming. The model is 
widely discussed and has two major limitations, no direct biological information and no minimal 
constraints to normal tissues. Linguistically, the purpose of radiotherapy may be stated as (a) 
delivering a desired tumor dose and zero dose outside the target volume; (b) delivering a high dose 
to the target volume and a low dose to the normal tissue. The statement (a) and (b) may be served 
as absolute linguistic prescription (ALP) and relative linguistic prescription (RLP), respectively. 
Although the ALP is ideal, it is clearly impossible to deliver due to the laws of nature [2], On the 
other hand, RLP clinically describes the strategy of radiation therapy. The words *high' and 'low' . 
used here are vague terms that are associated with the limitation of making precise definition. The 
complexity of treatment planning optimization is evident from the need to formulate some kinds of 
clinical goals to be optimized since there is no unique treatment plan which is clinically feasible 
and fulfills die two conflicting objectives: maximizing dose in the target volume while minimizing 
dose in normal tissues. 

Recently, several researchers [6]'[8] have come to pay attention to the analysis of 
\mcertainties in radiation treatment plaiming optimization. The tolerance of normal tissues has 
been discussed in [6] and [12]. Spirou et al [16] developed an inverse planning algorithm with soft 
constraints. The method allows acceptable doses of maximum and minimum as well as 
dose- volume constraints to the tissues of interest. 

The search for the optimal beams can usually be interpreted to be an optimization problem. Thus, 
the searching problem is converted to find the extremum of a given objective function. Several 
methods and algorithms have been investigated for inverse planning. Some examples are 
simulated annealing,^^ iterative approaches,^ as well as filtered back-projection and Fourier 
transformatioiL^'^*^ '^^ Although these methods are very promising, there are some aspects that can 
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be further improved upon. The simulated amiealing method may require long computation time 
due to the nature of random search. Most iterative approaches are parameter-dependent. The 
convergence and the quality of convergence may be affected by these parameters which are often 
determined by try-and-error. The filtered back-projection and direct Fourier transformation may 
have limitations on dose prescription and kernel selection. Inverse treatment planning is still at its 
early stage and many important aspects require be to further improved. 

In this research proposal, a new inverse planning method is introduced. A fuzzy approach 
is proposed to optimize the prescription of normal tissue. The presented method is based on the 
theory of fuzzy sets [9], and attempts to sufficiently use uncertain information under the tolerance. 
This new method contains two types of optimizations: intensity-modulated beam optimization and 
normal tissue prescription optimization. The former employs the fast-monotonic descent (FMD) 
technique. In this technique, a new iteration method is being developed in which the update 
scheme is analytically determined to avoid defected convergence. 

Execution of IMRT conformal plan using a dedicated delivery system requires accurate 
patient positioning. If patient is no correctly positioned, conformal radiation beams may be 
delivered to normal tissues rather than the planned target Therefore, patient mis-positioning may 
limit the applicability of dose escalation that is the key for MRT. 

A conventional radiation field is documented by use of a portal film in two-dimensional 
version. Information included in this image may not sufficient for IMRT procedxire, because the 
leaf position is not stationary during treatment for each field. Most quality assurance procedures 
for IMRT are performed in phantom. It is therefore, important to find a way to verify both 
anatomically and dosimetrically for IMRT treatment. It has been noted that monitoring actual 
dose delivered in IMRT using megavoltage computed tomography (MVCT) and portal imaging 
taken together with transit dosimetric method grows in great importance [Webb]. At present, a 
rapid and cost-effective method of verifying conformal IMRT radiotherapy based on limited 
number of fields is currently unavailable in clinical practice. 

In this proposal, a combined method will be developed to perform three-dimensional 
verification of patient setup and to document dose distribution treated using limited number of 
static IMRT fields. In this method, a megavoltage CT reconstruction technique will be developed 
based on Multilevel Scheme Algebraic Reconstruction Technique (MLS-ART) (Ref ?) using a 
megavoltage x-ray imaging device. By combining the transmitted treatment beams with the 
regular CT imaging projection beams, both patient geometry at treatment position and actual dose 
distribution can be reconstructed. The geometry and dose can be compared to the patient setup 
position and prescribed dose and used to correct subsequent beam placement or dose delivery 
accordingly. 

Portal CT and portal dose reconstruction is a novel verification technique in radiation 
therapy (especially in IMRT) with several advantages. It is online and allows direct verification of 
MRT for both patient position and dose delivery. Moreover, mega-voltage CT-based technology 
may potentially replace conventional patient simulation which uses kilo-voltage simulator or 
diagnostic x-ray CT and mega-voltage CT-based images may potentially be used for treatment 
planning. 
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C. PRELIMINARY STUDIES: 

C.l. Fast-convergence optimization algorithm 

The optimization of intensity-modulated beams (IMBs) consists of two main tasks: 
modeling (selection of objective function) and optimization (method of minimizing objective 
function). In this context, modeling means that the construction of a model in which our 
knowledge (physical, biological and clinical) about the irradiated structure's response to radiation 
is expressed by an objective function. The task of optimization is to develop a method by which 
one can obtain the optimal solution of minimizing the objective function. 

For multilateral optimization of radiation treatment planning, how to improve computation 
efficiency is an important topic [13], [17] and [18]. We are developing a general iterative method 
for intensity-modulated beam optimization. In this method, an optimal step-length, the key 
parameter in the update scheme for iteration, and an optimal solution to the problem of negative 
intensity are analytically derived. Therefore, the convergence to global minimum is not only 
guaranteed but also fast and monotonic descent. This method is called fast-monotonic descent 
(FMD) method, that can provide an optimal solution to the intensity-modulated beams either when 
the intensity valixe is greater than zero or when a negative solution is encountered. 

Let X = (^1,^2, ...,x^) be an intensity vector; Xn is the nHx component of intensity vector x . 
For each dose point (ij\k) , let P^^ represent the prescribed dose, and denote the calculated 
dose 

where >4^^ is a non-negative constant coefficient that can be directly calculated. The weight 
Wyi^>0 is used to indicate the importance of matching prescription and calculation. A qiiadratic 
objective function is therefore defined by 

/w=SZZ%*(^^-^.*)'- (2) 

i J k 

In the case of an optimization problem having an objective function of Eq.(2), the minimum cost 
problem is that of finding an achnissible intensity vector such that objective fimction /is niinimized. 
This constrained optimization problem can be written as 

(3a) 
(3b) 



minimizelfix)] 
subject to x^tH Vn . 
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Now consider an unisynchronous updating scheme used in iteration method. For an arbitrary 
evolution time /, when /->/ + !, 



x„il) + 6x„ if n=m 

(4) 

jc„ (/) otherwise 



and 

/(x(/))^/(x(/ + l)), 

where m is one of (1, 2,...,/2,...^. The updating scheme (4) says that, for each evolution time /, 
only one variable is adjusted. If each of variables is adjusted one time, then it is called one cycle, . 

Based on the theory of classical minimum, the necessary and suflBcient condition of 
descent for /is that tiie iterative rule satisfies: for each n 

^.=-lM. (5) 

where A„ is a small positive number and called step-length. Note that the iteration sequence 
generated by Eqns (4) and (5) is not guaranteed to converge to the minimum of/ This convergence 
is always dependent upon the choice of . Adequate selection of this parameter is critical for 
the success of iteration method. Generally, die choice of is a craft that is problem-specific. 

We found, for a quadratic fiinction, this parameter can be analytically derived and / will converge 
rapidly and monotonically to the minimimi with the following condition: 

^^ = 0 V„. (6) 
Parameter A„ can then be derived from the condition listed above. 

K •. (11) 

I J k . . 

With these two conditions (Eqns (5) and (6)), / descends rapidly to the global minimum if for each 
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N 



^«(0+Z5^,(/)+c„ 



^„(/ + l) = 



(7) 



Otherwise; 



where / denotes the /-th iteration, 



2 S S ^sr* ^mjjkA,yk 




i J * 



' J k 



and 



Z Z Z ^ff*^fli.(/Vt^ife 



< J k 



i J k 



The FMD algorithm can be summarized as follows: 

1) Fix the maximum number of iterations i, weights {w^j) , and termination criterion e >0. 

2) Initialize x(0) - (xi(0),X2(0),...,x^(0)) , and x^ > 0 for each n. 

3) For/=l,2,...,/:; 

a. Update intensity vector using Eq.(4). 

b. Compute = maxj|x„(/ + 1) - x„(/)|l . 




c. IF Ej<s stop; ELSE next /. 
4) Compute dose distribution using Eq.(l), 

Note that /is a constrained quadratic objective function. A set of values Jt^i ,^2 ^—y^N 
satisfies the non-negative constraints expressed by Eqn (3b) is called an admissible vector, and the 
admissible vector that minimizes the objective function is called the optimal admissible vector. An 
optimal admissible vector may fail to exist for two reasons. There are no admissible vectors (i.e., 
the given constrains are incompatible) or there is no minimum (i.e., there exists a direction in N 
space where one or more of the variables can be taken to negative infinity while still satisfying the 
constraints). Fortunately, neither of them is satisfied in the problem of intensity-modulated beam 
optimization. First, it is clear that, the sets in Eqn (3b) are convex, and the intersection consists of 
many points. Therefore, the non-negative constraints in Eqn (3b) are compatible. The second 
reason is also false since the intensity variables are non-negative. 

C.2. Fuzzy function for vague prescription 
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There is one important parameter {w^^} in Eq.(2) that has not been addressed in the 

previous section. Typically, we.know the prescribed dose for the target volimie and the upper limit 
for the sensitive organ. The prescription for the normal tissue is usually not given. Therefore, the 
optimization result varies with the prescription selected for the normal tissue. An intuitive strategy 
for finding the optimal normal tissue prescription would be to compare values of objective 
function calculated by using different prescribed doses and then to choose the minimum. In this 
way, w„ , the weight for the normal tissue, is a function of , the prescribed dose for the normal 
tissue. Here, the subscript n represents a point (1,7, K) inside the normal tissue. The difficulty of 
using this strategy is how to formulate the relationship between weight w„ and prescribed dose 
. Generally, all one knows is a plausible relationship between them: is the least when p^ 
approaches to zero and w„ is the greatest when p^ approaches to the upper limit. In this paper a 
dynamic weight function is used to express this fuzzy relationship. An optimal prescription dose 
for normal tissue is then determined by a validity function. 

C.2.1 Dynamic weighting 

Consider a quadratic objective function as shown in Eq. 2 with fiizzy weight as follows 



Pn\ if Uj>keQ„ 



and 



w,, if iJykeQ, 
w„ if iJykeCl^ 
w,. if hj\k^n„ 



p^ , p^ and p„ denote the prescribed doses for the target volume, the sensitive organ and the 
normal tissue, respectively. , and Q„ represent regions of these three corresponding 
structures, w^.^ e[0,l] is called fuzzy weight function that is used to emphasize the importance of 
matching the prescribed dose and the calculated dose for the point (/, 7, k) . Instead of fixing {i^^^^} 
and {w.jj^} in the hard inverse plarming (HDP), p„ is defined as a variable and w„ . is represented by 
a function of p^ in fuzzy inverse planning (FIP). Also, it is assumed that 



and 



I ' if Pn>Pn 
g{p„), otherwise 



(9) 
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where P, represents the prescribed dose in the target volume, /J is the tolerance dose in the 
sensitive organ, and is ttie tolerance dose in the normal tissue. g{p„) is a continuous function 
that increases with p„ . g(i>J = l when p„ is equal to the tolerance dose Here g(pj is 
called fuzzy weight function. 

Regarding the function g(') , one has only some vague knowledge that can be stated by the 
following two fiizzy rules: 1) the closer p„ is to , the closer is to one ( w„ =1 means the most 
important); 2) the closer p„ is to zero, the closer is to zero ( w„ =0 means the least important). 
Our purpose here is to use fuzzy technology to express this vague knowledge and to achieve an 
optimal solution. The form of fUzzy weight function may be obtained from planner's experience. 
As will be seen below, however, it is effective to use the following function: 

8iP,) = iff. 0 (10) 

where a: is a positive constant that controls the pattern of g(0 . These functions are shown in Fig.l. 
Obviously, any function with a K value of equal to or greater than 1 can be selected to express the 
mathematical meaning of the following linguistic prescription. For the normal tissue, the closer 
the prescribed dose p^ is to the tolerance dose, the greater the importance of the lose (i.e., the 
difference between the tolerance dose and the prescribed dose) is. 

C.2.2, Validity for inverse treatment planning 

Fuzzy inverse planning allows many feasible solutions to occur for a specific clinic 
problem. Selection of a specific treatment plan is determined by evaluating planning validity. Note 
that Eq. (8) cannot be served as a validity function since flie objective of radiation therapy 
optimization for non-target volume cannot be expressed by a quadratic function. To measure the 
validity of radiation treatment, we introduce a validity function, say v({P^J;x) , which is written 

as 

v({/^,};x)= Z IPifk-D^kl^ S D,j,. (11) 

Where, the first term represents the degree of dose uniformity for the target volume, and the 
second term represents the grade of protection of the non-target volume. For prescription validity, 
we flius presume that the ideal dose distribution can be achieved by minimizing v( {/Jy^} ; x) under 

the tolerance of normal tissues: 

minimize[v{{P^j^};x)]. (12) 

With the introduction of this validity function, the fuzzy inverse planning (FIP) algorithm can be 
summarized as follows. 

The FIP Algorithm: 
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FIPl . Fix the xnaximum number of trial dose prescriptions T for normal tissue, the 
maximum number of iterations i, and Ihe termination criterion e >0. Choose 
fiizzy weight function gi;) . 

FIP2. Initialize x(0) = (xi(0),X2(0),...,x^(0)), and > 0 for each n. 

FIP3.Forf=l,2,...,7 

1) Given (usually increases with r on the interval (0,i^]). 

2) Calculate w„ using Eq.(9). 

3) For/=l,2,..,L 

a. Update intensity vector using Eq.(7). 

b. Compute = max||x„(/ + 1) - x„(/)| . 

c. IF £^<f go to 4); ELSE next/, 

4) Calculate V using Eq.(ll). 

5) IF v,^, > V, stop; ELSE next t. 

6) Compute dose distribution using Eq.(l). 

C3. Evaluation of fuzzy inverse planning method using phantom cases 

The FIP is evaluated by two artificial examples. Dose-volume histograms (DVH) of the 
target volume (TV) and the sensitive organs (SO) are used as a primary tool for presenting and 
conq>aring dose distributions. 

C3.1. Casel 

This is a simulated cylindroid object aiid its central slice is illustrated in Fig.2. The 
geometry of this slice is similar to a CT head axial cut with two sensitive organs (analog to eyes) 
which are very close to the target volume. The prescription was given as follows: 100 dose units to 
the target volume, 20 dose units to the sensitive organs, and upper limit of 60 dose units to the 
normal tissue. Seven fan beams as shown in Fig.2 were unifomily arranged between 0-2 t: . Based 
on the primary-only model [10], the dose at depth is estimated by means of the percent depth dose 
data that were measured firom a field size of 4x4cm with 6MV photon beams: However, beam 
divergence was included. 

In order to show the convergence and fastness of FMD method, the FMD algorithm was 
run using three different values of the step-length: Jl„=O.Ol,yl^ = X^^^ and /l„ =0.001 
(;i=l,2,. . Here X^p, is the optimal step-length. Figure 3 shows the differences of convergence 
behavior between Xn^^opt and =0.001, and between A„ and ;i„=0.01 («=1,2,...,A0- 
This result indicates that the effectiveness of iteration methods is dependent upon the choice of 
step-length. The optimal step-length X^p^ derived here, however, provides an optimal 

performance in both the speed of convergence anid the quality of convergence. Figure 4 shows that 
FMD can provide a satisfactory result after 1 0 cycles. 

C Effect of normal tissue dose 
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In this study, a validity functioa is introduced to judge the optimal normal tissue 
prescription. Variation of validity function v versus prescribed normal tissue dose is plotted in 
Fig,5. The data in Fig.5 indicated that p„ =25 dose units appears to be the optimal prescription for 
normal tissue. Table 1 shows the fuzzy inverse planning (FEP) performance as a function of the 
normal tissue prescription dose with K=5 and Z=100. Here =0 means no normal tissue is 
considered in the FMD optimization algorithm, i.e., = 0 and = = 1 . Although statistic 
indices for the target volume in the case of p„ =0 are better tiian others, the average dose of 52.6 
dose units and the standard deviation of 33 .8 dose units for the normal tissue far exceed the upper 
limit of 60 dose units. The data listed in Tab. 1 shows that the optimal balance between objectives 
of high target dose and low normal tissue dose is achieved when p„ =25 dose units. Figure 6 shows 
corresponding dose-volume histograms for the target volume (Fig. 6(a)), the sensitive organs (Fig, 
6(b)) and the normal tissue (Fig.6(c)). The improvement of performance is evident with the 
optimization of normal tissue dose prescription. 

C.3.1.2. Comparison of FIP and HIP methods 

The case described above, with a normal tissue dose prescription of 25 units, can be 
considered an optimal result for the FIP algorithm: In this section the result obtained by the FEP 
method is compared to that obtained by the hard inverse planning (HIP) method. HIP means that 
only FMD algorithm, as discussed in Section n, is applied. For the HIP method, two extreme 
prescriptions are selected: (a) no normal tissue is considered in the optimization algorithm, i.e., 

w„ = 0 and = = 1 ; (b) the prescribed normal tissue dose is fixed, i.e., p„ = = 60 and 

w = w = = 1 . 

Dose-volume histograms are calculated and illustrated in Fig. 7(a) for the target volume 
obtained using FEP, HlP(a), and HIP(b), respectively. Corresponding dose-volume histograms for 
the sensitive organs and die normal tissue are illustrated in Fig. 7(b) and Fig. 7(c). Table 2 provides 
the relevant statistical parameters. It has been shown that the overall results obtained by optimizing 
prescription of normal tissue dose (FIP method) are better than those obtained by HlP(a) and 
HIP(b). 

C.3-1.3. Effect of parameter K 

As discussed in Section IE, validity function g( ) is considered to be adequate if K>1, 
The effect of the parameter K in Eq.(lO) on the performance of FIP method is evaluated by using 
values of K = 1,2, . ..,5 . Dose-volume histograms obtained using different K values for the target • 
volume (Fig. 8(a)), the sensitive organs (Fig, 8(b)), and the normal tissue (Fig.8(c)) are calculated 
and compared. Results indicated that the uniformity of dose for the target volume is improved as K 
increases. At the same time, the control of dose for the sensitive organ is stronger as K increases. 
However, the control of dose for the normal tissue is weaker as K increases. Therefore, the result of 
K = 5 is more desirable than those of others are. If K is greater than 5, the performance of normal 
tissue would be less desirable despite of dose improvement in other structures. 

C.3.2. Case 2 
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The central slice of the phantom geometry in Example 2 is illustrated in Fig.9. Similar to 
Example 1, seven fan beams are arranged at 0, 2n 11, An 11, 6 nil, 8 nil, 10 nil, and 12 nil, 
respectively. The prescribed doses were set: 100 dose units for the target volume, 20 dose units for 
the sensitive organs and the upper limit dose of 60 dose units for the normal tissue. The other 
parameters were Zr=100 and K=5. 

Figure 10 shows dose-volume histograms obtained by using FIP, HlP(a) and HIF(b) 
methods for the example 2. Figure 10(a) conesponds to the target volume, Fig.lO(b) to the 
sensitive organs, and Fig. 10(c) to the normal tissue. Table 3 indicates the relevant statistical 
parameters: means, standard deviations for the three structures in Example 2. The performance 
patterns of the FIP algorithm in Example 2 is consistent widi the results obtained in Example 1. 
However, the optimal normal tissue value here is equal to 30. Note that in this example one 
may also obtain a desirable result without considering normal tissue, i.e., FIP has similar result as 
HIP(a). 

C.4. M eavoltage CT image using limited number of projections 

We investigated using a fluorescent/CCD-based EPED, coupled with a novel Multilevel 
Scheme Algebraic Reconstruction Technique (MLS-ART), for a feasibility study of portal CT 
reconstruction (Ying 1990, Wong 1990, Yin 1994, Zhu 1995). We used an EPID, set it to work at 
the linear dynamic range and collimated 6 MV photons from a linear accelerator to a slit beam of 1 
cm wide and 25 cm long. We performed scans under a total of 200 MUs for several phantoms in 
which we varied the number of projections and the MUs per projection. The reconstructed images 
demonstrated that using the new MLS-ART technique, megavoltage portal CT with a total of 200 
MUs can achieve a contrast detectibility of -2.5% for an object of size 5mm x 5mm and a spatial 
resolution of 2.5 cm. 

We reported our investigation using a CsI(Tl) transparent scintillator x-ray detector 
together with the multi-level scheme algebraic reconstruction technique (MLS-ART) for 
megavoltage computed tomography (CT) reconstructions. The reconstructed CT images may be 
useful for positional verification in radiotherapy. The CsI(Tl) imaging system consists of a 
scintillator screen coupled to a liquid-nitrogen-cooled slow-scan CCD-TV camera. This system 
provides better contrast resolution than the standard electronic portal imaging system (EPID), 
which is especially useful given the low number of projections we are aiming at. The geometry of 
the imaging system has also been optimized to achieve high spatial resolution (I Ip/mm) m spite of 
the thickness of the screen. We present the images reconstructed using a pediatric head phantom 
using a total of 99 projections, and a combined phantom using 50 projections. Image 
reconstruction was carried out using the MLS-ART technique. We also present the CT images 
obtained using the back projection technique for comparison purposes. The detectibility? 

In addition, we also investigated the use of the kinestatic charge detector (KCD) combined 
with the multi-level scheme algebraic reconstruction technique (MLS-ART) for x-ray computer 
tomography (CT) reconstruction. The KCD offers excellent detective quantum efficiency and 
contrast resolution. These characteristics are especially helpful for applications in which a limited 
number of projections are used. In addition, the MLS-ART algorithm offers better contrast 
resolution than does the conventional convolution backprojection (CBP) technique when the 
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number of projections is limited. Here we present images of a Rando head phantom that was 
reconstructed by i^ing the KCD and MLS-ART. We also present, for comparison, the images 
reconstructed by using the CBP technique. The combination of MLS-ART and the KCD yielded 
satisfactory images after just one or two iterations. 

The advantages of MLS-ART applied to conformal radiotherapy are following: 

a. The MLS-ART outperforms the conventional CBP technique for low contrast detection given 
a limited number of projections and it is especially useful for megavoltage CT reconstruction 
since in radiotherapy we cannot rotate the linear accelerator gantry to acquire a large number of 
projections in a reasonably short of time. Contrast detectibility is strongly dose-dependent, and 
for some situations in x-ray imaging, high contrast resolution is not as important as the ability 
to provide excellent image contrast (Yaffe and Rowlands 1997). Such is the case with 
megavoltage CT imaging for radiation treatment verification. The high-energy x-ray photons 
experience inherently low attenuation in tissues. In addition, attenuation of radiation by tissues 
in this energy range is mainly due to Compton scattering that depends on electron density but 
not the atomic number. The two factors combined resulted in poor differentiation between 
various tissues (Johns and Cunningham 1983). Further, the detective quantum efficiency 
(DQE) of current megavoltage imaging devices is at least one order lower than those of 
detector for diagnostic x-ray CT. Therefore, megavoltage portal CT requires an efficeint 
reconstruction technique like the MLS-ART, especially the one that is optimal for situations of 
low-contrast detectibility. Better contrast detectibility also helps for more accurate dose 
reconstruction since spatial resolution imposed on dose is even more relaxed. 

b. MLS-ART can be used for CT reconstruction using the radiation treatment beams in addition 
to the regular CT projection beams. Such a 2-step reconstruction will produce much better 
reconstruction accuracy than simply using the treatment beams themselves because the latter is 
a case of incomplete data (although even for this case, MLS-ART itself works better than CBP.) 
In this way, the patient position can be directly and continuously monitored and even corrected 
during the treatment. 

c. Doing conformal radiotherapy using intensity modulated beams and portal CT is complicated 
by the tumor irregularity. Depending on the target shapes and sparing of critical organs, select 
treatment beam orientations to be orthogonal or close to orthogonal are important. These 
orientations must yield small geometrical correlations (less dose overlap) and most 
complementary dose distribution information. If possible, we could select the beam 
orientations following the MLS ordering, or in combination with the methods used by 
[Gokhale, Soderstrom, Bortfeld] . 

d. MLS-ART might be possible for dose reconstruction. It might be more accurate than the 
analytical method given a limited nxraiber of beams because for any analytical dose integration 
method, there is an implicit assumption that an infinite number of projections were used. But 
the integration method would fail if the angles between beams are large unless special 
techniques like arc therapy or tomotherapy are used. 
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C. 5- Reconstruction of IMRT beams for dose distribution 

D. EXPERIMENTAL DESIGN AND METHODS; 
D.l. Inverse treatment planning algorithm for IMRT: 

Aim 1.1: to develop a fast monotonic descent method for optimization of intensity 
modulated beams 

Aim 1 .2: to develop a fuzzy function to represent vague prescription in inverse treatment 
planning method 

Time table: 

Year one: Modeling of fuzzy function for vague prescription. Test of its efifectiveness xising 
three phantom cases: brain, pelvic, and head and neck. 

Year two: Implementation of inverse planning method into clinical practice. 

Year three: Test of inverse planning method with real clinical cases, 

D.l.l. Research method: 

First, it should been pointed out that the definition of ftizzy weight function as shown in Eq. 
(10) is not a unique form. Different functions may be used to achieve different goals. We tested 
Gaussian function instead of Eq. (10) and found that the result obtained using Eq.(lO) is better than 
that obtained using a Gaussian function. Second, in this study two loss functions are introduced. 
One is the objective function as shown in Eq. (8) and the other is the validity function as shown in 
Eq. (1 1). The former is used to optimize beam intensity. The latter is employed to evaluate the 
prescription of normal tissue. Note that the objective fimction as shown in Eq. (8) cannot replace 
the validity function as shown in Eq. (1 1). For example, it is equal important in terras of loss value 
when the calculated dose in normal tissue is either 10 below or 10 above the prescribed dose. 
Therefore, Eq.(8) does not completely express the objectives of radiotherapy, mid a validity 
function as shown m Eq.(l 1) is necessary. In addition, it is clear that Eq.(l 1) is an unbiased 
measure function since the loss for the non-target volume is calculated from zero. 

In this study we described a fiizzy inverse planning (FIP) method for solving the problem 
of uncertain prescription optimization in radiation therapy. This study concerned only with the 
optimization of normal tissue prescription. The dose prescription in the sensitive organs is fixed. 
Typically, the upper limit dose for the sensitive organ is less than that for the normal tissue. It is 
difficult to control the calculated dose less than the upper limit for the sensitive organ (except for 
those cases in which the sensitive organs are far from the target volume). Typically, the mean dose 
in the sensitive organs is greater than the upper limit dose. Note that, we assume that the 
importance of matching the calculated dose and the prescribed dose for the target volume is equal 
to that for the sensitive organ, i.e., = w, = 1 . Clinically, it means that the importance of 
protecting the sensitive organ is the same as that of controlling the tumor. However, different 
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weighting factors may be chosen by radiation oncologists for a specific clinical case to fulfill a 
special objective. 

A fuzzy inverse treatment-planning algorithm has been developed. This method provides 
an alternative to soft optimization for treatment planning. The main advantages have two folds, (a) 
The developed FMD has the fastest convergence speed in the stage of optimizing the beam 
intensity and the algorithm is simple to xise in which no parameter is problem-specific. And (b) the 
FIP technique can use uncertain information in inverse treatment planning to achieve the optimal 
balance between the objectives of matching the calculated dose and the prescribed dose for the 
target volume and minimizing the dose in normal tissue. The presented technique optimizes not 
only beam intensity distribution but also normal tissue prescription. The performance of the new 
algorithm has been compared to that of the hard inverse planning methods for two treatment 
geometries. The calculation time is less than 2 minutes on PC machine (333mHz, 64MB RAM) for 
10 slices with a matrix size of 256x256. At the present, it is difficult to compare between different 
approaches due to difference in test cases, dose calculation and other factors. It would be useful to 
unify the test geometry and to compare different methods in the future. 

Dl.1.1. Objective function 

Inverse planning method involves two key components: objective function to define the 
goal for the optimization, and optimization method to find the optimal solution for a given 
objective function. 

There are two categories of objective function: physical and biological. Since there are not 
enough clinical data to support most of biological models for TCP/NTCP, this study will continue 
to use an objective function based on physical models. The constraints for physical objective 
function involve dose uniformity and dose-volume requirements. We will use the quadratic model 
as described in the preliminary studies for establishing optimization objection. An analytical study 
indicated that some biological models could be approximated as a quadratic format (Ref.?). 
Therefore, the optimization algorithm developed here may be also applicable to resolve objective 
functions based on biological model. 

The physical meaning of wyk, Pyit, and Dqk are as described in the preliminary studies. 
D.1J.2. Optimization method 

The new iteration method, fast-monotonic descent (FMD) method as described in the 
preliminary studies has shown to be an very efficient approach. Therefore, we will continue to use 
this method for the optimization of intensity-modulated beams. 

Compared to some existing iteration techniques, there are several unique characteristics of 
FMD technique. (1) The key parameter, step-length, used in update scheme is analytically 
calculated so that no trial-and-error is involved. Choice of update scheme is critical for fast 
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convergence and optimal results. Inappropriate selection of the step-length may lead to poor 
convergence or even none convergence.- (2) Fast and monotonically convergence guarantees the 
global minimum of the optimization algorithm. (3) The problem of negative beam intensities is 
effectively eliminated. (4) The algorithm is simple to understand and implement for clinical 
applications. 

D.1.1.3. Fuzzy representation of vague prescription 

The concept for a logistic plan is to deliver full dose to the target region while keeping the 
dose below the maximiun tolerance for normal tissues. The proposed quadratic model in the above 
section is not sufficient to address the upper limits for normal tissue prescription. A fuzzy function 
will be introduced to represent vague normal tissue prescription. 

The theory of fuzzy set is a mathematical tool used to represent uncertain or partial 
knowledge. In inverse treatment planning, one only knows the upper limits for normal structures 
but is not certain what is the optimal prescription, especially the prescription for normal tissue and 
critical organs. Therefore, this category of problem can be represented by the theory of fiizzy set. 
We will develop an appropriate fuzzy function to represent vague or uncertain knowledge. 

The objective function will be divided into three terms. The first term relates to the target 
volume that is expected to receive uniform prescription dose. The second term relates to the 
critical organs that are sensitive to radiation damage and a tolerance dose will be set. The third 
term relates to the normal tissues except critical organs in which dose is expected to be as low as 
possible. To achieve these goals, the weight factors in the objective function will be redefined as 
following: 

Wijk in target volume Wt = tt, where tt is a constant and is used to indicate the importance of 
matching target dose. Typically, tt=l. 

Wijk in critical organs Wc = cc * h(pc,PcX where cc is a constant and is used to indicate the 
importance of matching prescription for each critical organ. 

Wijk in critical organs Wn = nn * g(pn,Pn), where nn is a constant and is used to indicate the 
importance of matching prescription for normal tissues. 

Here both g and h are two fuzzy functions. The fundamental of constructing a fiizzy 
function is to find proper weighting factors in the objective functions. Linguistically, the closer 
the prescribed dose is to the tolerance dose, the greater the importance of the lose (i.e., the 
difference between the tolerance dose and the prescribed dose) is. Mathematically, it may be 
described by a following function: 

8iPn>ify> ^^Pn^K (10) 

h(Dc) = (^)"' OSDcSPc (10) 

Pc 
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where Kis?i positive constant that controls the patterns of g(0 and /jf.) . These functions are 
shown in Fig. 1 . Obviously, any function with a K value of equal to or greater than 1 can be selected 
to express the mathenaatical meaning of the following linguistic prescription. Both exponential, 
and Gaussian function will be examined for this purpose. 

D.1.13. Penalty of optimization method 

Fuz2y inverse planning allows many feasible solutions to occur for a specific clinic 
problem. Selection of a specific treatment plan is determined by evaluating planning validity. Note 
that Eq, (8) cannot be served as a validity function since the objective of radiation therapy 
optimization for non-target volume cannot be expressed by a quadratic function. For example, 
when the calculated dose is greater than (but closer to) the prescribed dose, the objective function 
will not able to penalize such a situation. To measure the validity of radiation treatment, we 
introduce a validity function, say v({P^.J;x) , which is written as 

(11) 

The importance of requiring the quadratic function is that it is proved that the global 
minimum is exist and is unique. Here tt, nn, and cc are used to indicate the importance of matching 
each term. These parameters may be determined by the planner based on the clinical needs for 
each individual patient. When an equal inq>ortance is reached, tt, nn, and cc are equal to 1 . If 
DVH is used to judge the results, tt, nn, and cc can be used changed to reach final plan. 

D.1.2. Evaluation of inverse planning method: 

D.1.2J. Phantom test 

We will test the developed inverse planning method with three cylindrical phantoms: brain, 
head and neck, and pelvis. Each geometry has both complicated target volume and critical organs 
around it. Clinical applications: based on tiie infonnation available here, we will try to use 
step-and-shot technique to deliver the radiation. 

Several measures will be calculated for quantitative evaluation of each test : 
Dose volume histogram 

Statistical data: mean, standard deviation, minimum, and maximum. 

Initially, we will use primary beam only for the dose calculation. Beam divergence will be 
considered and inverse square law will also be applied. The primary beam is acquired from TMR 
data for 4x4cm field size of 6 MV photon beam. 

D.L2.2. Patient case test 
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D.1.2.2. IMRT experiment 

Clinical implementation will based on step-and-shot approach. 

a. patient CT image 

b. input to Pinnacle 3-D planning systme 

c. contour target and critical structure and external edges 

d. export contours and CT images to inverse planning algorithm 

e. generate IMRT intensities for each beam 

f. segment each beam for step-and-shot 

g. import segmented field to Pinnacle 3-D system 

h. calculate MUs for each segment 

D.2. Three-dimensional reconstruction of dose distribution 

Aim 2. 1 to develop a megavoltage CT reconstruction device 

Aim 2.2 to develop a dose reconstruction based on limited IMRT radiation beams. 

d.2. 1 Research method 

d.2.1.1 
d.2.1.2 

d.2.2 Research design 

d.2.2.l 
±222 



Time table: 
first year: 
second year: 
third year: 

With the patient at the treatment position, we will use the same projection data for the 
geometry reconstruction to estimate the true dose delivered to the patient. We will develop a scheme 
for 3-D dose verification, which requires overlaying the reconstructed patient geometry at treatment 
with the distribution of the delivered dose. This will serve as a verification tool to the initial treatment 
planning. Any major discrepancy may warrant modification of the patient setup and/or the treatment 
strategy. 
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In the current megavoltage CT imaging, we used uniform beam for both the calibration runs and 
the projection measurements. For the IMRT, the x-ray beam from each projection will be modulated 
in intensity, i.e., non-uniform. Therefore, we need to measure the IMRT beams before the patient is 
placed in the treatment room to get the entrance intensity distribution. An alternative way is to 
down-load this distribution from the MRT delivery files; however, this option is less direct than 
measurement Without these entrance beam intensity, we would be unable to decide whether any exit 
intensity change is due to the entrance intensity change or different attenuation within the patient 
geometry. Image reconstructions using the intensity modulated beam can be tested initially using 
simply compensator or wedge. 

We plan to use different patient dose calculation metiiods based on the measured transmission 
x-rays (most of the detected x-rays are primary components for that there is a '-SO cm air gap between 
the patient exit surface and the detector. The scatter fraction for a 20x20 field size, 17 cm thick water 
and 30 cm air-gap is 10% [Jaflfray et al., 1994]). The first two methods are based on the primary 
photon fluences at the point of dose calculatibii. We call them the convolution-superposition and the 
superposition-convolution method. In the first method, the x-ray fluence at the detector surface is 
ray-traced back to inside the patient's geometry. We convolve the fluence inside the patient with an 
appropriate energy deposition kemel (dose spread array) to obtain the dose distribution. By 
superimposing all the distributions over the reconstructed patient geometry, we can obtain the total 
dose distribution. In the second method, first, the x-ray fluences of all the beams at the detector 
surface are ray traced back and superimposed together to get the total primary x-ray fluence 
distribution inside the patient Then we convolve this total fluence distribution with a rotation dose 
spread kemel to get the total dose distribution. In both methods, normalization (calibration) is needed. 
The third method is based on the primary photon distribution attenuated inside the patient. The 
overall primary attenuation distribution in the patient, which is different from the total primary 
fluence ray traced back from the detected x-rays, may also be reconstructed, using the similar 
methods for the emission tomographic reconstructions such as PET and SPECT. (For each beam the 
attenuation profile is obtained by subtracting the penetrated primary from the entrance beam-) Then 
we convolve this overall attenuation distribution with a special rotation dose kemel (which is 
calculated based on the photon numbers rather than the photon fluence) to get the dose distribution. 
MLS-ART may be applied for such a photon attenuation reconstruction (with some modification) 
based on the experiences of Herman [1993], but not the conventional CBP technique due to the 
limited beam numbers. Further, the inhomogeneity corrections can be directly calculated based on the 
geometric reconstructions. Compared to the first two methods, the 3rd method may be more direct, 
accurate and convenient. It is also easier for the intensity modulated beams. With faster and growing 
computation technology including hardware specifically designed for MLS-ART and FFT, we may 
achieve faster and more accurate on-line verification. 

To be more accurate, we need to calculate the primary fluence on the detector's surface by 
deconvolving the projection data (measured during treatment session, therefore no additional dose to 
the patient) using the "EPID kemel" (the point spread function of EPID). 

We can also use some portal dosimetry methods to model the exit dose distribution and to 
compare it to the calculation results. We can also use it to model the absorbed dose inside the 
treatment volume (actually the dose along the beam path) based on flie treatment geometry as 
reconstmcted using the NCLS-ART. We will compare our results to some oflier dose modeling and 
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verification methods such as the portal dose imaging (PDI) technique [Wong et al., 1990], the 
superposition/convolution method [McNutt et al., 1996a & b], and the inverse filtered (convolution) 
back-projection method of Hobnes and Mackie [1994]. 

The other advantage of using the CsI(Tl) detector for dose modeling and verification throughout 
the treatment volume is that compare to commercial EPID which oveiresponse to low energy x-rays 
for dosimetry studies, fliis detector is more tissue equivalent If we use the a-Si detector, a more 
active way to reduce over-response is to use organic scintillators (low-Z plastic materials) on the 
detector so that detector response will be more tissue equivalent. One way is to use a low-Z screen 
with a buildup phantom such as the solid water. Hien the photons undergo interactions in the buildup 
material. The secondary electrons are mainly absorbed in the organic screen, and the visible photons 
are emitted toward the a-Si photodiode sensor. The merit of using such an organic screen is that the 
dose deposition by electrons inside the tissue can be exactly modeled by using the tissue-equivalent 
organic material. This new screen needs to be fabricated by a medical imaging company because that 
the currently available screen has poor surface smoothness. 

For clinical application of megavoltage portal CT, improving the accuracy of reconstruction 
rests on more efficient detectors and optimized reconstruction algorithms to most effectively use 
the available dose. In this proposal, three interrelated specific aims will be perused. 

(1) Adapt an efficient x-ray detector to carry out this project. We have 3 options: a) Using Varian 
Portal Vision Electronic Portal Imaging Device (EPDD), which tested to be the best 
commercial EPED system, b) Using the amorphous silicon system. Varian would lease us a 
piece of the detector for research purpose for 3 months, c) Collaborate with Herber Zeman's 
group in the University of Tennessee, Memphis to use their CsI(Tl) CCD system which was 
specifically designed for megavoltage imaging. This CsI(Tl) system is one of the best 
megavoltage imaging system providing both good contrast and spatial resolution. 

(2) Image reconstruction.. We will use the MLS-ART technique for this specific application, in 
which we first use a limited number of cone beam projections (dosage close to that in 
diagnostic CT) to get megavoltage CT reconstruction for patient geometry. Then we will use 
the treatment beams to further modify the reconstruction. This two-step reconstruction has 
three important purposes: 1 . the second stage will locally improve the CT image quality inside 
the tumor; 2. The second stage would also let us know the placement of treatment beams inside 
the patient geometry obtained by the first step. 3. The treatment (dose covered) area can be 
visualized from the final CT images. MLS-ART can easily perform this 2-step reconstruction. 
However, it is impossible for the conventional convolution backprojection (CBP) technique. 

(3) Dose reconstruction and verification. From the portal image taken at the treatment portal, we 
could obtain the portal (transit) dosimetry and convert this portal dose information to photon 
fluence. We then trace this fluence back to the patient and get the fluence inside the target We 
then use the convolution-superposition or superposition-convolution or other methods to 
calculate the dose inside the patient use the fluence. Such a 3D dose distribution can be 
overiaid onto the patient geometry to verify the treatment plan. Any major discrepancy 
between the prescribed and actual dose can be corrected by modification of the treatment 
setup. 
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In first step, we will test Aim (2) (image reconstruction) based on the work Huaiqun Guan already 
did (2 projects) on megavoltage CT reconstruction using the Phihps SRI-100 system and the 
CsI-CCD system. We will test Aim (3) using computer simulation. 

B. HUMAN SUBJECTS: 

No human subject will be used in this study. 

C. VERTEBRATE ANIMALS; 

No animals will be used in this study. 

D. CONSULTANTS: 

No consultant will be required. 

E. CONSORTTOM ARRANGEMENTS: 

No consortium will be arranged. 
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1. Describe the questions to be answered by the research. Include a title and an 
abstract as well as a general background statement similar to the introduction 
section of a scientific paper on the proposed work. This should include 
references, especially those of the principal investigator. 

Title: Artificial Intelligence-guided Inverse Treatment Planning 

Abstract; An artificial intelligence (AI) method, fuzzy logic, is applied to optimize 
parameters in the inverse treatment planning for intensity-modulated radiation 
dierapy (IMRT). With the capability of fuzzy inference, the parameter modification 
of the objective fimction is guided by physician's treatment intention and experience. 
For the different parameters involving inverse planning, the corresponding fuzzy 
inference systems (FISs) are developed in order to accomplish the treatment 
requirement. With the fimction of fuzzy inference, the efficiency and quality of 
inverse planning will be substantially improved. 

Introduction: 

The conventional inverse treatment planning involves optimizing an intensity 
spectrum by an objective function in order to obtain an ideal dose distribution. 
Sometimes, the dose distribution calculated by the optimized intensity spectrum might 
not be the one vrt expected. It is often resulted by an improper selection of parameters 
selected in the objective function. As the desired dose distribution can only be achieved 
by a refmed combination of parameters that are usually unknown, it was inevitable for a 
plarmer to try different combinations of parameters in order to achieve a desired dose 
distribution. At present, this process is done by a trial-and-error approach.^^'^*'' It is not 
only time-consuming but also difficult to find a desired combination of parameters. 

A number of methods have been proposed in recent years in order to tackle this 
problem.^"*'^*^^"^^^ Starkschall^*^^ proposed an approach that removed the necessity of 
defining a "best" treatment plan, and incorporated the dose-volume constraints into a 
system to search for a feasible plan that could satisfy the constraints. If no calculated 
doses satisfy the treatment goal, the planner will provide a guide about how the dose- 
volume constraints may be modified to achieve a feasible result. This approach is only 
applied to the conventional three-dimensional (3D) treatment planning. Wu and Mohan^"*^ 
developed an optimization system, which employed both dose- and dose-volume-based 
objective functions. In their system, the optimal plan is selected by calculating the cost of 
the objective function, or "plan score" (the lower the score, the better the plan). Xing et 

[13,14] presented a method .which employed a second stage evaluation function to 
compute the differences between the calculated and the ideal dose volume histograms. 
Based on the results of the evaluation function, the weighting factors in the objective 
function are adjusted. This procedure is to minimize both the objective and evaluation 
functions in a round-robin manner. Later, further improvement is achieved by using a 



statistical measure called preference function, which is constructed based on the 
empirical judgment. The problem about the selection of weighting factor still exists 
because it translated to the problem of how to specify the parameters in the evaluation or 
preference function. A similar method was also proposed by Wu et al.^^^' using a genetic 
algorithm to optimize the weighting factors and beam weights in the conventional 3D 
treatment planning. Li and Yin^^^ introduced fuzzy logic into the inverse planning system 
to adjust the weighting factors for normal tissue. The result was promising. However, 
optimizing the parameters for the target and critical organ were not included in the 
system. Also, the weighting factors initialized by the fuzzy functions still need to be 
modified by the trial-and-error approach. 

One of the major advantages of IMRT is to spare dose of the critical organ 
adjacent to the target. When a strict compromise between target dose and critical dose 
was required, fine-tuning of parameters in the objective function is important and 
necessary. As the effect of parameters on the output of the inverse planning is uncertain, 
it is hard to say which combination of parameters is best for the current case. Generally, 
this decision is made according the priority of organs involved and is individualized. 
More concerns will pay to the organ with high priority. As the trial-and-error approach 
can be regarded as a way using empirical knowledge to guide the parameter modification, 
it can be substitute by existed AI approaches, such as expert system and fuzzy inference 
system. The main purpose of this research is to use fuzzy interference system to define 
parameters in the objective function based on treatment intention instead of trial-and- 
error approach. 
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L INTRODUCTION 

Conventional inverse treatment planning involves optimizing the beam intensity 
spectrum by using an objective function in order to obtain an ideal dose distribution. 
However, the dose distribution calculated by the optimized intensity spectrum might not 
be the one we expected. The choices of weighting factors in the objective function have 
an important impact on the calculated doses. Usually, this kind of tasks is done by a trial- 
and-error approach based on planner's experience. ^"^^ It is not only time-constmiing but 
is also difficult to achieve the expected dose distribution among different organs 
involved. 

A number of methods have been proposed in recent years in order to tackle this 
problem."*'^*^ Starkschall^ ' proposed an approach that removed the necessity of 
defining a "best" treatment plan, and incorporated the dose-volume constraints into a 
system to search for a feasible plan that could satisfy the constraints. If no calculated 
doses satisfy the treatment goal, the planner will provide a guide about how the dose- 
volume constraints may be modified to achieve a feasible solution. This approach was 
only applied to the conventional three-dimensional (3D) treatment planning. Wu and 
Mohan"* developed an optimization system, which employed both dose- and dose- 
volume-based objective functions. In their system, the optimal plan was selected by 
calculating the cost of the objective function, or "plan score" (the lower the score, the 
better the plan). Xing et alP'* presented a method which employed a second stage 
evaluation function to compute the differences between the calculated and the ideal dose 
volume histograms. Based on the results of the evaluation function, the weighting factors 



2 



in the objective function were adjusted. This procedure was done to minimize both the 
objective and evaluation functions in a round-robin manner. An improvement was 
proposed using a statistical measure called preference function, that was constructed 
based on the empirical judgment.'^ Further investigation is expected to resolve the 
problem of how to specify the parameters in the preference function. A similar method 
was also reported by Wu et a/." using a genetic algorithm to optimize the weighting 
factors and beam weights in the conventional 3-D treatment planning. Li and Yin' 
introduced fiizzy logic into the inverse planning system to adjust the weighting factors for 
normal tissue. The result was promising. However, optimizing the weighting factors for 
target and critical organs were not included in the system. 

As inverse treatment planning plays a key role in IMRT, effectively selecting the 
parameters for a given objective function becomes a challenging problem. The 
increasing appUcations of IMRT has prompted the need for the development of an 
intelligent technique to guide the optimization of parameters used in inverse planning. 
As the fuzzy inference system (FIS) is constructed to perform the inference as a human 
does, it is suitable to solve problems involving parameter optimization in inverse 
planning." In this study, we developed a fuzzy inference system to guide the 
optimization of weighting factors in inverse planning. The principle of the fiozzy 
inference system is introduced in Sec. II. The results witii two representative cases (one 
simulated case and one clinical case) are presented in Sec. III. In Sec. IV, we discuss the 
strength and capability of this approach. 



n. MATERIALS AND METHODS 

For a given dose prescription, conventional inverse treatment planning consists of 
two steps: (1) finding the suitable weighting factors for involved organs and (2) 
optimizing the intensity spectrum based on the given weighting factors. As there are a 
large number of choices for weighting factors, finding the desired ones for a given 
objective function is difficult. In this study, we developed a fuzzy inference system (FIS) 
to optimize the values of weighting factors and compute the intensity spectrum by the 
iterative gradient technique. The involved organs in this system are categorized as the 
target volume (TV), the critical organs (CO), and the normal tissue (NT). 

A. The principle of the fuzzy inference system 

The flow chart of FIS is illustrated in Fig, 1 . It consistes of three main modules, i.e., 
the Fuzzifier, the Inference Engine (consisted of fuzzy rules) and the Defuzzifier. For 
each variable input to the fuzzy inference system, a number of fiizzy sets are defined with 
appropriate membership functions. These membership functions are labeled with 
linguistic tags frequently used by humans (such as "High" dose). During the process of 
fuzzification (corresponding to the module of Fuzzyifier), the single input value is 
compared to the membership functions defined for that input variable. If the membership 
functions have a non-zero output, it will take effect on the final result of the FIS. 
Generally, the fuzzifier calculates the response of rules for the input values, and the 
inference engine modifies the consequent of rules in response to the input values and the 
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defiizzifier generates final output based on the result of inference engine. 

The inputs of this system are defined as the characteristic doses [€„ , C^o , C„ ], 
which consist of the mean dose {Mean, , i =rK. CO. NT) and its standard deviation 
{STD„ i ^TV, CO, NT). For the target volume, C^= Mean^-STD„ . For the critical 
organs and normal tissue, Cco= Afea«co + 5rZ)co and C^= Mean„+STD„. The 
outputs of FIS [ APP'jv , AWco , ] defined as the adjustment of the weighting 
factors for each involved organs. For each input variable, two fiizzy sets, "High (H)" and 
"Low (L)", are defined with membership fianctions [f" (x) , /^^ (jc) , / =TV. CO, NT\, as 
shown in Figs. 2a-2c. For each output variable, three fuzzy sets, "Increase (I)", "Not 
change (N)", and "Decrease (D)", are defmed with membership functions 
[ g' W . zf W . &f (*) » ' ^ • *^8et volume, these three membership 

fimctions are shown in Fig. 2d. Similar membership fimctions are defmed for critical 
organ and normal tissues for the same adjustment strategy and therefore not shown. 

Based on the input and output variables defined above, fiizzy rules are established 
for the fuzzy inference engine. As shown in Appendix I, eight rules are employed in this 
system. In each rule, the if-part of rule is called antecedent and the then-part of rule is 
called consequent. Two of them (Rule 5 and Rule 8) are used to demonstrate the 
procedure of fiizzy inference as shown in Fig. 3. Note that the input (output) variables 
are labeled using the bold fonts and their corresponding linguistic tags are labeled using 
the italic fonts in each rule. According to the linguistic tags, the corresponding 



5 



membership functions for the input fuzzification are specified as shown in Step 1 of Fig. 
3. Such as the input variable C^y , the membership function /^^ is specified in Rule 5 by. 
the linguistic tag "High". For each rule, the outputs of the fuzzification are [ , D^o , 
] and [D^ , D^o , D^j- ], respectively. Based on these outputs of fuzzification, the 
degree of support ( Auppoct) ^^^^ achieved by a logic operator "Min", such as 
A'upport=Min(Z)^,, D'co. Dlr)md Nippon = Min (i)^,, Z)*^, Z)]^), as shown in Step 2 

of Fig. 3. Note that the degree of support represents the applicability of rule's antecedent 
for given mputs. Based on the degree of supports, the fuzzy inference is performed by a 
standard implication method, which is accomplished by a logic operator "Min". For 
example in Step 3, the membership function (x) in Rule 5 is modified as 

grf' W Alpport I Stv (^) )• Th^ modified membership functions became 
[g^;-^ gcf'i^)^ W] and [g^^^ W, gl:^ix), for Rule 5 andRule 

8, respectively- As there are two sets of modified membership functions obtained, it is 
necessary to combine them to produce a smgle one. In Step 4, they were aggregated into 
one set by a logic operator "Max", i.e., [grvix) , g^o W, g^(x) ]= [MaxCg^^-^^ (x) , 
g^^W],Max(g^;f^(x), gi:'^ix)\M?.xig'/^ix), g'^{x))l The aggregated 
functions represent the combined consequent from all the rules. Finally, the aggregated 
functions are defuzzified to a single value by the centroid method in Step 5. The x- 
coordinate of the centroid (represented by sign "©") for each aggregated function was the 
final output, the adjustment amount of weighting factors. 
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B. The fuzzy logic guided inverse planning algorithm 



The flow chart of the FLGIP system is schematically illustrated in Fig. 4. First, the 
dose prescription and weighting factors are set to their initial values. Then, an iterative 
gradient algorithm is used to calculate the intensity spectrum x. hi our study, the 
objective function is defined as follows: 



/(x) = Z Z Z (Pifi -^iikf^ 



(1) 



/ ; k 



Where d^f^ = calculated dose for each voxel, A„ y^ is the relative dose 

net 

coefficient, or dose per unit intensity of pencil beam, py^^ is the dose prescription and 



w^jf^ is the weighting factor defined as follows: 



Pco^ if (iJ^)^^co> 



Wjy, if iij\k)eCl^ 
Wco. if (ij\k)enco^ 

^NT^ if iiJ>^)^^irr 



Cl^ , and denotes the target volume, the critical organ volimie, and the normal 

tissue volume, respectively. The minimization of the objective function imder the 

constraint of x„ ^ 0 can be written as a problem of 

min{/(x)} 

subject to x„ > 0,Vw. 

Eq. 2 can be solved by the fast-monotonic-descent (FMD) method developed by Li and 
Ym^ which is an optimized iterative gradient technique for the quadratic function. Based 
on the optimized intensity spectrum, the characteristic doses are calculated and then input 



to the FIS. Using fuzzy inference, the adjustment amounts of weighting factors 
[MV^,6iWco' APT/yj.] are obtained. Then, the weighing factors for the next iteration are 
modified as follows 

As the weighting factors affect the output of inverse planning by their relative values 
rather than the absolute values, they are re-normalized to [0, 1] by the following formula: 

w;in . 1) = ^^"^^^ — / . {TVMNT} . (4) 

This updating procedure repeats until the following convergence criterion (5) is satisfied: 
j[C^ (n + 1) - C^in)Y ^ [Ccoin + 1) - C^q jn)]' + [C^jn ^ 1) - CMf < t (5) 

where T is a small threshold number, such as 0.01. 



IIL RESULTS 



The performance of FLGIP system was examined using two cases (one simulated 
and one clinical). Dose-volume histograms (DVHs), plus the variation of characteristic 
doses and weighting factors versus the iteration number, are used as the primary tools to 
evaluate the performance of this system. Pencil beams of 6 MV were used. For 
simplicity, the primary-only dose at depth is used in the calculation. The initial 
weighting factors [^^k.^co^^/vtI are set to [1,1,1] (after re-normalization using formula 
(4), they became [0.58, 0.58, 0.58]) and the convergence constant T was set to 0.01 . 
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A. The Simulated Case 

The central slice of this case is illustrated in Fig. 5a. The layout on this slice 
simulates the spinal cord with a target volume surrounding it. Seven treatment beams are 
uniformly arranged between 360 degrees. This configuration is typical in spinal 
radiosurgery using IMRT. 

The FLGIP system was tested using four sets of different dose prescriptions: [100%, 
20%, 50%], [100%, 30%, 50%], [100%, 40%, 50%], [100%, 50%, 50%]. Fig. 6 shows 
the variation of (a) characteristic doses and (b) weighting factors versus the iteration 
number for dose prescription [100%, 30%, 50%]. The results indicate that for the target 
volume and critical organ, the characteristic doses monotonically converge to the 
prescribed doses (the normal tissue dose also converges, but in a much less rate due to its 
large volume.) The results shown in Table I demonstrate that the high target dose and 
low critical organ dose are achieved simultaneously and both meet the prescribed ones. 
Their corresponding DVHs for (a) the target volume, (b) the critical organ, and (c) the 
normal tissue are shown in Fig. 7. Note that, the final results also depend on the provided 
dose prescriptions. For each set of dose prescriptions, their corresponding isodose 
distributions are shown in Fig. 8. 

The effect of initial weighting factors on the final characteristic doses was examined 
by using eight sets of initial values with the same dose prescriptions [100%, 30%, 50%]. 
The characteristic doses for each set converged v/ithin 50 iterations. The final results and 
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the standard deviations are shown in Table II. The results indicate that the achieved 
characteristic doses by different sets of initial weighting factors are comparable. We 
averaged the final eight sets of weighting factors. The mean weighting factors and their 
standard deviations are 0, 139 i 0. 11 3 for the target volume, 0.985 ± 0.025 for the critical 
organ, and 0.004 ± 0.003 for the normal tissue. It implies that there could be a limited 
space for the selection of optimal weighting factors. 

B. The Clinical Case 

The central slice for this clinical case is illustrated in Fig. 5b. Eleven beams are 
arranged at 0^ 33 ^ 66°, 90^ 120^ 150", 210°, 240°, 270°, 300°, 330°. This 
configuration represents a complicated IMRT case. The dose prescription is set to 
[100%, 30%, 50%]. The variations of the characteristic doses and weighting factors 
versus the iteration number are shown in Fig. 9a and Fig. 9b, respectively. The 
characteristic dose CrKmonotonically converges to its prescribed dose 100% while the 
characteristic doses Ceo and C^vr monotonically converges to the doses below their 
prescribed values,'30% and 50% respectively. The DVHs of the calculated doses for 
different organs at three iterations 5, 10, 15 are shown in Fig. 10. The results indicate that 
the gap between the DVHs of target volume (Fig. 10a) and critical organs (Fig. 10c) 
increases with increased iteration number. The substantial improvements of isodose 
distributions around the critical organ COi (the one closet to target volume) in different 
iterations can be easily identified from Fig. 1 1 . 
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IV. DISCUSSIONS 

A fuzzy inference system was developed to automatically modify the weighting 
factors in inverse treatment planning in order to achieve the dose distributions best . 
matching ttie treatment requirements. The fundamental inference mechanism is 
demonstrated by a mini system consisted of two rules as shown in Fig. 3. Among the 
eight rules, Rule 5 plays the primary role to drive the system toward the convergence 
while Rule 8 (plus the other six rules) drives the inputs toward its prescribed ones. For 
example, when the inputs for critical organ and normal tissue are much higher than their 
prescribed doses, the output of FIS will mainly be determined by the adjustment of Rule 
8. Once flie inputs approach their prescribed ones to better match the antecedent of Rule 
5 (usually after several iterations), the consequent of Rule 5 will take more effect on the 
output of FIS and drive the system towards convergence. Note that, the other six rules in 
Appendix I are used to process different scenarios of mismatching between characteristic 
doses and prescribed doses of different organs. 

The details of adjustment process are shown in Fig. 6 and Fig. 9. At the first several 
iterations, the weighting factor for the target volume decreases and the weighting factor 
for the critical organs increases quickly. After few iterations, as the characteristic doses 
approach their prescribed ones, the adjustments of weighting factors gradually reduce. It 
is noteworthy that the characteristic doses for the target and critical organs in the last 
iteration satisfy their dose prescriptions. For the normal tissue, however, the final 
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characteristic dose is appreciably lower than its prescribed one due to its large volume. 
Although some rules seemingly take less effect on or are seldom used in these two cases, 
these rules are necessary for the more complicated cases. In addition, the results shown 
in Fig. 7 indicate that using different dose prescriptions could result in different dose 
distributions. Potentially, the fuzzy inference technique may also be used to optimize 
other parameters in inverse planning such as the beam orientation, the dose prescription 
etc. Investigations are still underway for these applications. 

As the configuration of FIS is flexible, it provides us a wide space to customize the 
configuration for different applications. In this system, the input characteristic doses are 
chosen as the mean dose combined with its standard deviation. For target, the lower than 
mean input dose will help the FIS to drive the target dose to be higher toward the 
prescribed one in the next iteration. Similarly, for critical organ and normal tissue, the 
higher than mean input dose will drive critical and normal tissue doses to be lower 
toward the prescribed ones in the next iteration. In this way, both high target dose and 
lower critical organ (and normal tissue) doses may be easier to achieve. For output 
variables, they are simply defined as the relative adjustment of the weighting factors, 
which are between -1 and 1 . For the selection of inference rules, it is primarily 
determined by the clinical experience. The general treatment mtention may be described 
as: If the target dose is low, its weighting factor should be increased. If the critical organ 
and normal tissue doses are high, their weighting factors should be increased. In our 
system, such treatment intention is expressed by eight rules (Appendix I), which is a 
complete combination of linguistic tags for three kinds of involved organs. This option 
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may not be the best bui can avoid any unpredicted input values. As for the selection of 
membership functions, the Gaussian function is adopted due to its simplicity and 
popularity for most of the engineering applications. In some circumstance, part of the 
Gaussian function is used, such as those shown in Fig. 2a-2c. 



V. CONCLUSION 



A fuzzy logic guided inverse planning system has been developed. This system 
provides an effective and efficient approach to optimize the parameters used in inverse 
planning. The main advantage of using FIS is that it can perform the sophisticated 
inference formerly done by trial-and-error approach. Relymg on the planner's experience 
and knowledge on how to compromise parameters among dififerent organs involved, the 
optimization of weighting factors can be easily accomplished by encoded rules. As 
demonstrated by the result of two cases, the fuzzy inference system could undertake the 
very complex task of parameter optimization in inverse planning. 
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FIGURE CAPTIONS 



Figure 1 Schematic illustration of the fuzzy inference system (FIS) used for modification 
of the weighting factors. 

Figure 2 Illustration of membership functions used in FIS. (a) The membership functions 
"High"(/^(;c))and"Low"(/^(x))definedforinputvariable C^. (b)The 
membership functions "High" (/coW) and "Low" (/^(;c)) defined for input 
variable Qo . (c) The membership functions "High" {f^{x) ) and "Low" 
if NT W ) defined for input variable . (d) The membership functions 
"Decrease" (g^(x) ), 'T^ot change" {g^{x)\ and "Increase" (g^^(x) ) defined 
for the output variable tiWjy . 

Figure 3 Demonstration of the inference procedure by two rules. Step 1 : Inputs 

fuzzification. Step 2: Degree of support. Step 3: Fuzzy inference (implication 
operation). Step 4: Aggregation operation. Step 5: Output defuzzificaiton. 

Figure 4 Flow chart of the fuzzy logic guided inverse treatment planning system. 

Figure 5 The central slices for (a) a simulated case and (b) a clinical case. TV, CO and 
NT represent the target volimie, the critical organ and the normal tissue, 
respectively. Arrows pointed to the target volume indicate the beam directions. 

Figure 6 The variations of (a) characteristic doses and (b) weighting factors versus the 
iteration number in the simulated case for the dose prescription [100%, 30%, 
50%]. Note that the initial weighting factors [1,1,1] were normalized to [0.58, 
0.58, 0.58] using formula (4). 
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Figure 7 Dose-volume histograms of the calculated doses in the simulated case for (a) the 
target volume, (b) the critical organ, and (c) the normal tissue for four sets of 
dose prescriptions. 

Figure 8 The dose distributions in the central slice of simulated case for four sets of dose 
prescriptions (a) [100%, 20%, 50%], (b) [100%, 30%, 50%], (c) [100%, 40%, 
50%], (d) [100%, 50%, 50%]. 

Figure 9 Variations of (a) characteristic doses and (b) weighting factors versus iteration 
number in the clinical case for the dose prescription [100%, 30%, 50%]. Note 
that the initial weighting factors [1,1,1] were normalized to [0.58, 0.58, 0.58] 
using formula (4). 

Figure 10 Dose-volume histograms of calculated dose distributions for five involved 
organs: (a) the target volume, (b) the normal tissue, (c) the critical organ 1, (d) 
the critical organ 2, and (e) the critical organ 3 at iteration 5, 10, and 15, 
respectively. 

Figure 1 1 The dose distributions in the clinical case at (a) Iteration 5, (b) Iteration 10, and 
(c)Iterationl5. 
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ABSTRACT 



An intelligent technique, fuzzy logic, was applied to optimize the weighting factors 
in the objective function of inverse treatment planning for intensity-modulated radiation 
therapy (IMRT). Based on this technique, the optimization of weighting factors is guided 
by the fuzzy rules while the intensity spectrum is optimized by a fast-monotonic-descent 
method. The resultant fuzzy logic guided inverse planning (FLGIP) system is capable of 
finding the optimal combination of weighting factors for different anatomical structures 
involved in treatment planning. This system was tested using one simulated (but 
clinically relevant) case and one clinical case. The results indicate that the optimal 
balance between the target volume and the critical organ dose was achieved by a refined 
combination of weighting factors. With the help of fuzzy inference, the efficiency and 
eflFectiveness of inverse planning for IMRT are substantially improved. 
Key words: fiizzy logic, optimization, inverse treatment planning, intensity-modulated 
radiation therapy. 
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